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Abstract 

We calculate quark-antiquark potentials using the relationship between the expectation 
value of the Wilson loop and the action of a probe string in the string dual. We review 
and categorize the possible forms of the dependence of the energy on the separation between 
the quarks. In particular, we examine the possibility of there being a minimum separation 
y—i for probe strings which do not penetrate close to the origin of the bulk space, and derive a 

condition which determines whether this is the case. We then apply these considerations to 
the flavoured resolved deformed conifold background of Gaillard el al [T]. We suggest that 
the unusual behaviour we observe in this solution is likely to be related to the IR singularity 
which is not present in the unflavoured case. 
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1 Introduction 



One of the most important observables in a field theory is the Wilson loop , given by 

W(C) = j r JrPe'i A , (1) 

for a closed curve C. In particular, the potential E of a quark-antiquark pair is related to the VEV 
of a Wilson loop. Indeed, for a separation L, take the loop C to be a rectangle of sides L, T, with 
T — > oo, then 

(W(C)) ~ e- £r (2) 

The Wilson loop is particularly significant in the context of gauge-string duality and the Mal- 
dacena conjecture [3J, because it is accessible from the string side of the correspondence. For a 
review see [3]. As proposed in [5J[B], the associated quantity in the dual string theory is the action 



of a string world-sheet which ends on C at the boundary of the AdS space ( Figure 1 1 . That is 



W(C)= f VFe~ s[F \ (3) 

JdF{C) 

where F describes the fields of the string theory, with boundary values dF(C). In the limit of 
strong coupling the result is that the Wilson loop corresponds to the area of a surface bounded 
by C, extending into the bulk and forming the world-sheet of a classical string. This means that 
in the strong coupling regime of the QFT 

(W(C)) ~ e" 5NG , (4) 

and referring to |Equation 2[ the energy of the quark-antiquark pair corresponds to the Nambu- 
Goto action of the string |H [7] 

E= S f. (5) 

This is divergent, and is renormalised by subtracting the (infinite) quark masses, given by the 
action of two rods from the ends of the string to the end of the space, as descibed in [3 HI |7] . This 
is discussed more rigourously in [5]. 




Figure 1: Schematic diagram of how the Wilson loop C (blue) relates to the world- sheet of a string 
extending into the bulk. The loop should be taken to extend an infinite distance in the t direction. 
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The aim is therefore to solve the equations of motion resulting from the action and so determine 
the shape of the string formed for a given L in a given background. In |section 2| we discuss the 
equations of motion for a string, following the derivation in 01 [7] , and describe a generalisation 
to D branes. We discuss the possible behaviour of the function E(L), and of the string shape 
in sections |3](4j We then demonstrate the results of the preceding sections with respect to some 



well- understood backgrounds ( section 5 1 , before applying them to the flavoured resolved deformed 



conifold [T] in section 6 



2 The action and equations of motion 
2.1 Action 

We will consider backgrounds of the form 

ds 2 = g^dx^dx" = -g tt dt 2 + g xx dx 2 + g pp dp 2 + g, ,dO'd(l> , (6) 

where gu, g xx and g pp are functions only of p. Here and throughout this paper we use the string 
frame. We will restrict our attention to p-dimensional objects which extend on time and one 
spatial Minkowski direction x, and probe the radial direction according to 

x = x(X 1 ), P= P (X 1 ), (7) 

where we use world- volume coordinates X a , < a < p. For p > 1, the object also extends in the 
internal space described by the coordinates (9\ We are interested only in the static case, so we 
can identify X° = t. 

The action for such an object, with tension To, is 

S = T / d" +1 Xe- a V-detG Q 0, (8) 

where 

G a f) = 9nv gXa (9) 

is the induced metric on the world volume, and a = 1 for a Dp brane, with a = otherwise. For 
the configuration described, the induced metric is 

^Luccd - -9tt(dX ) 2 + (g xx x' 2 + g ppP ' 2 ) (dX 1 ) 2 + G^ l) dX a dX\ (10) 



Gab — 9ij Q Xa Q Xb i (11) 



,(p-i)_ d0< d9i 

*ab ~ 9*3 Q Xa g X , 

where x' — dx/dX 1 . Writing the time interval as T, the action (Jsj) is then 

S = T tJ dPX e- Q * y/g tt ( 9xx x' 2 + g pp p' 2 ) detG^ . (12) 

This can be written in a form corresponding to a 1-dimensional 'effective string' with tension 
Tstr(p)- Defining f(p) 2 = g tt g xx and g(p) 2 = gu9 PP , 



where 



S = T IdX 1 (fJ x ' 2 + ^p' 2 ) T str , (13) 



T Btr (p) = T e-«* J d^X y/detG^-V. (14) 
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This has a simple interpretation: The tension T str is as expected the energy density on the effective 
string, while the factor 



dX 1 f 




(15) 



is the length element on a string embedded according to ^ in the geometry ([6]). The action 
( 13 1 is therefore an obvious generalisation of the case considered in jH [7J to a string with p- 
dependent tension. In fact ( 13 ) can also be obtained from the action used in [7J section II. A] by 
the replacements 



/ -> T stI f, 



9 -> T stI g. 



(16) 




Figure 2: The generic shape of the effective string (solid line), 
solution (251, discussed in subsection 2.3 



The dashed line shows the 'free' 



2.2 Rescaled radial coordinate 

The following discussion will be considerably simplified by the introduction of a rescaled coordinate 
R defined by 



dp / 



Then the metric ^ becomes 



f 2 

ds 2 = -g tt dt 2 + J — (dx 2 + dR 2 ) + .gydfl W , (18) 

gu 



and the action ( 13 1 is 

S = T J dX 1 Vx' 2 + R' 2 T cff , (19) 

where T c r = fT str . Notice we can interpret this action as being that of a string in any metric of 
the form 

ds 2 eS = -Jttdt 2 + lxx {dx 2 +dR 2 ), (20) 
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provided we give the string a tension equal to T e g/ y/'yttlxx- In particular, interpreting T c g itself 
as the string tension results in Minkowski space. The behaviour of the effective string is then 
described completely by the function T c r. This is in contrast to the description in terms of the 
original metric in equations (|6| and ( 18 ), for which the natural interpretation involved a geometric 



factor (15 1 as well as an independent tension T str . 

Although useful in the general discussion, the integration involved in obtaining R(p) means 
that this coordinate system will be difficult to apply to any specific case except for extremely 
simple backgrounds. 

2.3 Equations of motion and separation 



The derivation of the equation of motion from (JT9J) is essentially the same as the calculation in 
[4, 7 . Imposing time independence, we get the single equation 

dR , dx 



C dx^ ± dx^- C2 - < 21 > 

Parametrising the effective string as 

X 1 = x, R = R{x), (22) 

we can integrate to obtain the shape of a string with a minimum radial coordinate at R = Rq. 
We impose C = T cff (i? ), and obtain 

, ™ ] 0<x<^, 

..■<*. = r- ^ vw ; r # )2 " " 2 (23) 

dR' , T M) ^ = , ±<x<L, 




where 

L(Ro) = 2 / dR 

is the separation of the endpoints of the effective string and i?oo = R(p — > oo). This is the same 



L(R ) = 2 (*~ dR-=^^== (24) 



as is obtained by modifying the result of [7] according to the prescription (16 1. 

In some cases we will find that T(Rq) — 0, where Rq is the minimum radial coordinate contained 
in the space. Then there is an additional solution to pTJ (with C = 0) which is not compatible with 



the parametrisation ( 22 ) . This corresponds to a string which drops vertically from the endpo ints 
and stretches horizontally along the 'bottom of the space', R = &o, as shown in Figure 2 A 
suitable parametrisation is 

f(o, i^o-^ 1 ), x x <o, 

(x,R) = I (X\ Ro), Q<X X <L, (25) 

{(L, Ro-L + X 1 ), X 1 > L. 

As we shall sec, gcncrically L(Rq) has inversion points, and together with the possibility of 



the extra solution (25 1 this means that Rq(L) can be multivalued. The different branches can 



be interpreted as corresponding to stable, metastable and unstable configurations for the effective 



string (see section 3 1 



2.4 Boundary conditions in the UV 

When we consider a fundamental string we must enforce Dirichlet boundary conditions at R — > 
i?oo, as described in [7J. This corresponds to the string ending on a D-brane at large R. Specifically, 
we require that 

^0 (26, 
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for R — > Roq. Referring to (21 1 this means that we need 

A. — — = °°- (27) 

Recalling that we have imposed C = T c ff(i?o) and that for a fundamental string T e g = f, this 
becomes 

hm = oo. (28) 

R^R^ f(R Q ) 2 

As this must hold for all Ro, we can simply require that 

lim f(R) = oc. (29) 

R—^Rod 

Although this condition is not required when the string is replaced by a D-brane, we will 
restrict our attention to those backgrounds in which (29 1 holds. 



3 Energy and stability 



As noted in 



section 1 



the energy of the quark-antiquark pair is simply given by S/T, which 
as would be expected corresponds to the tension integrated along the string. This is in general 
infinite, and so we renormalise by subtracting the action of two vertical rods extending from Rq 



to infinity [51 !H HI E] • For the smooth solution (23), which we can parametrise as R(x), this gives 



E{R ) = 2j^ dRdl 



dxV 

— r ff - 2 
dR) 



R„ 



Using the equation of motion (21 1, this can be written 



E(R ) = 2 dR T cfi (R) 
JRo 



T oS (R) 



^T eS (R)2 - T cff (i? ) 2 



dR T cff . 



(30) 



2 1 dRT eS (R). (31) 

'Ro 



Given a form for the function L(R ) it is simple to obtain the qualitative behaviour of E(L) 



without evaluating (31 1 
force is 



Generalising the result obtained in [71 [S] using ( 16 ) we have that the 



dE 
dL 



= T eS (R Q ). 



(32) 



Given the reasonable assumption T ff(i?o) is continuous and positive this implies that the extrcma 
of E(Ro) correspond to those of L(Ro), and to cusps in E(L). 

The possible presence of extrema in the function L(Rq) raises the question of which branches 
of the solution represent physical configurations of the string for a given L. It was shown in |10j 
that these extrema also correspond to the boundaries between stable and unstable configurations 
(ignoring any regions of instability due to fluctuations in the angular directions l ). Although this 
implies that only one side of the extremum describes a physical solution (as opposed to one being 
stable and the other metastable), it is not clear how to identify the physical branch. However, we 



can make progress if we assume that T e g(Ro) is always increasing with Rq. Referring to (32 1, this 
is equivalent to the statement that 



dR 



[#(L(iZo))] >0. 



(33) 



In terms of the function E{L), we see that E"{L) changes sign at each cusp (Figure 3 1, with 
E"(L) > for the upper (higher E) branch and E"(L) < for the lower. We can relate this to 



G 



E 




Figure 3: A generic cusp in E(L). The concavity condition leads us to expect that the upper 
(dotted) branch is unstable. The arrows show the direction of increasing Rq. 



the concavity condition discussed in [11] . namely that for a physical interaction between quarks 
we must have E"(L) < 0. We therefore expect that the upper branch at each cusp is unstable. 

Probably the simplest form of behaviour occurs when T e g(i?o) 7^ 0. By the argument of 
[T2l [3] this means that E(L) becomes linear at large L, corresponding to confinement. If L(Rq) 
is decreasing for all Rq, we obtain the qualitatively simple behaviour exhibited, for example, by 
the Klebanov-Strassler [12] and wrapped D5 brane [H] models. This is depicted schematically in 
|Figure 4| 




Figure 4: The qualitative behaviour of (a) L(Rq) and (b) E(L) in a simple confining case, such as 
Klebanov- Strassler . 



More complicated behaviour is exhibited in the models with massive dynamical flavours |15| 
[TSJ [TTJ , and also the walking D5 background discussed in detail i n (19117]. There L(R ) has two 
local extrema, leading to two cusps in E(L), as shown in Figure 5 We still have T(Ro) ^ 0, and 
confinement is again seen for large L. In (7 an analogy with a van der Waals gas was proposed. 
This again suggests that the upper branch at the cusps, corresponding to L'(Rq) > 0, should be 
identified as unstable. 

When T c s(Ro) — we also obtain the second solution (251. As we renormalise by subtracting 
the action of two vertical rods, and there is no contribution to the energy from the part of the 
string with R = R , this solution has E = independent of L. This is the stable solution for 
sufficiently large L, so at large separations the endpoints of the string behave like free particles. 

As pointed out in [TU], this is analogous to the case of a soap film stretched between two circular 
rings. In fact, if T c ff(i?) cx R, corresponding to a string in Rindler space, the analogy becomes 
exact [2D], as the action (19 1 is then identical to that of the soap film. The 'free' solution (251 
corresponds in the case of the soap film to a disconnected configuration, with the film forming a 
disc over each of the rings. 
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Additionally, when i?o = Ro the integrand in (24 1 vanishes for all R ^ Rq. This means that 
unless the lower limit of the integral gives a non-zero contribution the separation given by the 



smooth solution (23 1 will go to zero as the string approaches the end of the space; L(Rq) — 0. As 
can be seen from Figure 6[ this can be considered a special case of the 'van der Waals' behaviour 
(Figure 5 1, in which the minimum in L(R ) moves to the origin. 



(a) 



(b) 



Figure 6: The qualitative behaviour of (a) L(R ) and (b) E(L) in the 'soap film' case. The red 
region now corresponds to the solution (|25[). 



If the lower limit of the integral in (24 1 diverges, the separation will diverge for i?o — > -Ro despite 
having T e ff(-Ro) = 0. This is the case for a string in AdSs x S 5 . As before, we can consider this a 
special case of the 'soap film' case, with the maximum in L(R ) moving to (R Q — Rqi L — oo), as 
shown in Figure 7 This results in a qualitatively Coulombic potential (exact in the case of AdS.). 
The 'free' solution is now presumably metastable for all L. 

Finally, it is possible for any of these forms of behaviour to be further modified by the devel- 
opment of an additional local maximum and minimum, as in the 'van der Waals' case, |Figure 5] 
This will again result in a pair of cusps. 



4 General results on the behaviour of L and E for large Rq 

In the discussion of |section 3| it was implicitly assumed that L — > for Rq — > R^ . This is not 
always the case. We will find cases in which instead we have L — > ^ 0. 
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Figure 7: The qualitative behaviour of (a) L(R ) and (b) E(L) in the 'Coulomb' case. 



Assuming that T(R ) — > oo for R a — > Roo, we can write the separation (24 1 as 

T 



L(R ) 



dT 



1 



T'{R) Ty/T*/T% - 1 ' 

where in this section we write T = T c s(R) and Tq = T ff(i?o)- Defining t = T/Tq, this is simply 



(34) 



/OO 



t 1 



in which the Rq dependence is entirely contained in the factor 

t = T(R(t,R )) 
t'(R) T'(R(t,Ro)y 



(35) 



(36) 



where the factors of Tq from the definition of t cancel. 

We are interested in the case when the separation is constant for large Rq. The integrand in 
( 35 1 must therefore be a function only of t, which is equivalent to requiring that ( 36 ) is a function 
only of t for large R. That is 



T 



T'(R) 



T 



for R^R„ 



(37) 



However, the left hand side is explicitly independent of Rq so this can only be satisfied if T(t) is 
a constant. We therefore have L(Rq — > R^) = constant if and only if 



T 



T'{R) 



constant 



for R^ R B 



In this case the integral (35) can be evaluated, giving 



Lao = lim L(Rq) = 7r lim 

-Ro-s-floo R^R, 



T 



T'(R) 



(38) 



(39) 



The above calculation of relied on the fact that the integral for L(Rq) covers only the 
range R < R < R^. The integral (31 ) for E(R Q ) covers the whole range R < R < i?^, and so 
an analogous calculation is not possible. When L(Rq — » R^) — > constant, it is however possible 
to find a condition which determines whether E(Rq — > Roc) — > constant. 

Referring to ( 39 1 , for large R we can write 



T 
T'(R) 



'-'OO 

7T 



A'(T(i?)), 



(40) 
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where A'(T(R)) -> for R -> i?^, so that by d35 



L(Ro) — L c 



dt A'(T(t,R )) 



The energy is given by (31 1, which becomes 



E(R ) 





^+A'(T) 




7T 



1 



)- 2 /> 


^+A'(T) 




7T 



where To = T c ff(i?o). After evaluating some integrals this can be written as 

T 



or T" 1 

T(i? ) = — — + 2A(f ) - 2A(T ) + 2 



dT A'(T) 



vt 2 - n 



i 



(41) 



(42) 



(43) 



The first two terms are constants, and the last two, which contain the Rq dependence, involve only 
the region Rq < R < i?^. Whether the energy approaches a constant is therefore determined by 
the large R behaviour of A'(T(R)). The condition is that E(Ro) — > constant for Ro — > R^ if and 
only if A'(T) vanishes at least as fast as 

A'(T)~^, e >0 (44) 

for T^oo. 

The generalisation of the discussion of |scction 3] to account for non-zero L x is simple. In 
general there is a region with L < in which the 'free' solution (25 1 is the stable one, as in 



Figure 8 If there is a minimum, so that L(R) approaches the asymptote from below, then there 



is an additional (presumably unstable) branch as in Figure 8 (c,d). 

However, in some cases L(R) is always increasing, as in Figure 8| (e,f). In this case the con- 
siderations discussed in |scction~3| suggest that the 'free' solution is the only stable one for all L. 



5 Application to specific cases 



The results obtained above are not immediately useful when we are interested in a specific back- 
ground, as opposed to general considerations. This is because, as was noted in |subscction 2.2| the 
integral 



R(p) = J" dp' 



f(p') 



(45) 



is either difficult or impossible to evaluate analytically in all but the simplest cases. The results 
of sections [2]|4] are easily generalised to allow the more general coordinate p, by the insertion of 
factors of 

dR _ g 
dp ~J' 

In particular, 

g(p) T cff ( Po ) 



(46) 



L(p ) =2 dp 

J Pn 



PO 
oo 



f(p) V T M 2 ~ T ofi(Po) 



r 



(47) 



E(p ) = 2 dp 

J Pa 



g(p)T c e(p) 
f(p) 



Tett(pa) 



yr cff ( P )2 - T cff ( Po )2 
Similarly, ( 39 ) is more conveniently written 

Tes(p) g(p) 



"d,»f. (48) 



7r lim 



(49) 
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5.1 The Klebanov-Strassler model 



To demonstrate these ideas, we will look at the case of a string in the Klebanov-Strassler back- 
ground |13| . This is a convenient, relatively simple, example of a theory which results in confining 
behaviour for large L. 

If we restrict our attention to 1-dimensional objects, we need only the non-compact part of the 
metric, 



ds 2 = h'^dxl 



1 



4/3 



h 1 ' 2 



dp 2 



where 



h(p) = a 



2 2/3 



x coth x — 1 , , „ .1/0 
dx — (sinh2x- 2x) 1/3 , 



sinh x 



K(p) 3 = 



sinh 2p — 2p 



2 sinh 3 p 

The functions appearing in the action (l3| are 

fipf 
9(P) 2 

while for a string T e g = /. 
For small p, we have 

h = h — h 2 p 2 + • • • , 



(50) 

(51) 
(52) 



£ 4/3 

6K 2 ' 



(53) 



.9 



j7 = .9o + 92P 2 + 



Then the contribution from the lower limit of the integral (47 1 is 

goh 



1 



I'u 



Vp 2 -p 2 o + --' 



so that L diverges logarithmically for p — > 0. Using (32), we find that for small po 

E'(L) -> T cff (0) = constant, 



(54) 



(55) 



(56) 



as expected for a confining theory. 

The results described in |scction 4| have limited application to this case, simply telling us that 
L does not approach a non-zero constant for large p. The relevant function ( 49 ) is for large p 



T e s(p) g{p) g{p) 



T'M Hp) f'( P ) 



Vpe 



-p/3 



(57) 



which is not constant. 

To obtain the behaviour over the full range of p it is necessary to integrate ( 47|48 1 numerically. 
The result is shown in |Figurc 9 The expected confining behaviour is seen for large L, and the 
form is qualitatively that of Figure 4| 



5.2 AdSVSchwarzschild x S 5 



An example which produces the 'soap film' behaviour (Figure 6 1 is the large mass limit of the 
ylrfS'5-Schwarzschild black hole, expected to describe finite-temperature JV = 4 SYM |21j . This 
was discussed in [2"2"] . 
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The metric is 



ds 2 



K 1 



dt 2 + dx 2 



n 2 



dr 2 +ft 2 dft 2 , 



so that g(r) = 1 and 



(58) 



(59) 



The horizon is at ?o = /i, and we will consider strings, for which T e ff = f. Then T c g(r) — > for 
r — > fo, as is necessary for the'soap film' behaviour. 

This background is in fact simple enough that we can obtain the rescaled coordinate R(r) 
exactly, although not in a form which is particularly useful. We have 



K 2 



dR 



and so we can define 

R(r) = K 2 

This results in R(?o) = and 



dr' 



/tFT i 



K 2 TC 1 



sjr' 4 - n 4 



1 1 _ 5 ji 4 

4' 2' 4' V 4 



(60) 



(61) 



Roo = R(r -4 oo) 



(|) ^ 2 
r(|) m 



To find the behaviour for large r we need the function 

T(r) g = TZ 2 Vr 4 - fi 4 
T'{r) / 2 r 3 



2r ' 



Referring to section 4 we sec that this will result in L(ro) — > for tq — > oo. 
We can write the separation in the form of ( 34 1 , as 



L(r ) = UT cS (r ) 



dT r 



eft 



Toff(ro) 



Ti 



ft 4 



-3/4 



1 



VT e 2 ff -T cff (r )2 



(62) 



(63) 



(64) 
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For small r, T e s(r) — > 0, and so the lower limit of the integral contributes 

K 3 1 



L(r ) = TZT cS (r ) 



TMro) ^ V T cS - T cff (r ) 



T c g ( r ) log T off (r ) 



It is actually possible to evaluate ( 64 1 exactly, resulting in 

'3 3 5 



2^r(f) 



ft 2 



L(ro) r (if (r 4 - ^i 4 ) 1 ^ 2Fl ^4' 4' 4' r 4 — p* 



(65) 



(66) 



Together with the function E(L) , obtained numerically, this is shown in Figure 10 The relationship 
with |Figure 6| is clear. 



(a) 



(b) 



0.8 
0.6 
0.4 
0.2 



To r °-2.5 




1.0 



Figure 10: Plots of (a) L(rg) and (b) E(L) for a string in the AdS'5-Schwarzschild x S 5 background 



(581, for n = TZ=l. 



6 The flavoured resolved deformed conifold 

Having discussed the necessary techniques and general results, we can turn to the main new 
material presented in this paper. The analysis of [23] describes a system of solutions related by a 
chain of dualities, together with a boost in eleven dimensions. More specifically, the 'unrotated' 
solution, obtained by setting the boost parameter (3 = 0, corresponds D5 branes wrapping the S 2 
of the resolved conifold |24j , and is a simpler limit of the solution with general /3 as introduced in 
|25j . This has additional D3 brane charges which are not present in the 'unrotated' case. 

Taking the limit f3 — > 00, we obtain the 'rotated' solution, which describes the baryonic branch 
of the Klebanov-Strassler theory. 

The metric is of the form 

ds 2 = h- 1 ' 2 dxl t!i + e 2 *h 1 ftd4. (67) 
The warp factor is related to the dilaton by 

h = e" 2 * - e~ 2 *~ tanh /3, (68) 

where is the asymptotic value of <I> for large p. 

These expressions also apply in the flavoured generalisation, described in [1 . However, the 
/o-dependence of the functions is different. Most significantly, the solutions are now singular in the 
IR: $ — > —00 for p -» 0. 

Turning to the UV, we find that to leading order Q(p) is unchanged by the addition of flavours. 



However, the form of (68 1 is such that in the 'rotated' case (/3 = 1) is such that h(p) is sensitive to 



14 



the sub-leading behaviour of The 'rotated' solution therefore has different UV asymptotics in 
the flavoured case. This is interpreted in [T] as resulting from smeared source D3 branes, uniformly 
distributed in p. These result from the action of the 'rotation' on the source D5 branes, in the 
same way as in the unflavoured case the 'rotated' solution has bulk D3 branes resulting from the 
colour D5 branes in the 'unrotated' solution. 

The aim here is to use the methods discussed in the previous sections to assess the physical 
significance of these changes with respect to the unflavoured solutions. The changes in the UV 
asymptotics can be isolated by considering strings which do not descend close to p = 0. The 
results derived in |section 4| will therefore apply. 

It appears more difficult to isolate the effects of the IR singularity, as the string always probes 
the large p region as well. However, in most cases we will find that the effective tension T e g 
vanishes for small p. As discussed in |section 3) this means that only the lower limit of the integral 



(24 1 contributes to L(p). The limiting behaviour of L(p ) for small po is therefore presumably 



insensitive to changes to the UV. 
6.1 The solutions 

We will now define the solutions of interest more concretely. The metric is of the form 
ds 2 = h-V 2 dxl s + e™h l ' 2 dsl, 

dsl = e 2k dp 2 + e 2q {iol + u%) + \e 2 $ [(w x + a^) 2 + (u> 2 - au 2 ) 2 ] + ^e 2k (Cj 3 + uj 3 ) 2 , (69) 
where 

uji = d9, Qi = cos ip d9 + sin tp sin 9 dip, 

lu 2 = sin 9 dip, ui 2 — — sin ip d9 + cos ip sin 9 dip, 

lu 3 = cos 9 dp, &3 = dip + cos 9 dip. 

The coefficient functions {<&, g, k, q, a} depend only on p, and as above 

h = e- 2<s> -e- 2 *~tanh^. (70) 

The coefficient functions were shown in |57] to be given by 



e 



2<I 



1 ~ *> = P*»kT-Q, a= (71) 



4Pcoshr-Q' ^' Pcoshr-Q' 

and the BPS equations reduce (after some choices for constants of integration) to 

sinh t 



sinh 2p 



Q = 2 ^ f (2/Q cosh r — 1), 



e 4(*-*o) 



(P 2 -Q 2 )e 2k sinh A T' 



e 2k = l(P' + N i ), (72) 



together with 



The solution discussed in [T] is given by the asymptotic behaviour of P{p), 



h lP +±f -plogp - £plog(- logp) + O ( p} °^ p] ) +0(p 3 log P ), 



P = 



ce 4p/3 + 9N, 



(2N c -N l ) 2 (p 2 -p+ 1 jl)~ 8 -^ e-^ + 0(pe^), p -4 oo, 



(74) 
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where the two arbitrary constants hi and c arc related in a non-trivial way. 

The full solution to ( 73 1 can then be found numerically, interpolating between the two regimes 
in (l74|- 




Figure 11: The generic behaviour of (a) P(p) and (b) Q(p), obtained numerically. The red 
and green dashed curves show the first few terms of the expansions (74) for small and large p 
respectively. 



6.2 The 'rotated' case 



The functions appearing in the action (13 1 are 



e- 2 * -e- 2 *~tanh/3' 



9(P) 



2 _ 2*+2fc 



(75) 



To obtain f(p—> oo) — > oo, as required by the boundary conditions for a string (subsection 2.4 1 we 
therefore require the limit (3 — > oo. This is the 'rotated' solution; the field theory limit discussed 
in p. 

6.2.1 Fundamental string 

In the case of a fundamental string we simply have T c n = f /2wa'. We first consider strings which 
descend deep into the space. For small p, the functions we need have the asymptotic behaviour 
[TJ Appendix B] 



o 4($-* ) 



2k 



27 

2Nf 



(-logp)" 



1 + 



log(-logp) 



(-logp) 



1 + 



- log p 
log(- log p) 
-logp 



Writing 27ra' = 1 this results in 



T eS (p) = Hp) 



1/4 



e*°(-logp)- 3 / 4 



1 + 



log(-logp) 

-log/3 



(76) 
(77) 

(78) 



In particular, T c ff(0) = 0, so we generically expect qualitatively 'soap film' behavior as in Figure 6 



or one of the modifications discussed subsequently. The 'free' solution (25), with po = 0, exists 

smool 

L(p ) = 2 I dl 



for all L. In the case of the the smooth solution (23 1 the separation is given, as in section 4 by 

1 



t 



t'{R) tVt 2 - 1 ' 



(79) 
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Generalising this to allow us to continue working with p rather than R, we can write 



L(po) = 2 



dt 



t g(p) i 



(80) 



where we have split the integral at an arbitrary point p — p* As T o ff(0) = 0, the limit po —> 
results in 



t(p*) = - , { -> oo. 



The separation is then 



L(p o ^0) = 2^ 



Tefi(Po) 

t g(p) 1 



dt 



t'(p) f(p) tV^T' 



(81) 



(82) 



This integral covers only the range < p < p* . As was arbitrary, we can take the limit — > 
and evaluate (82 1 exactly using the small- p asymptotic expression for tg/t'f. Then, using ( 76|77 l, 



t'(p) Kp) f'(p) 



-p(-logp) 3/2 



1 + 



log(- log p) 

-l0g/9 



We also have 



so 



f(p) /logp \ 3/4 



f(Po) V lo gP 



and 



L(Q) 



4 /2iV f 



lim p (-logpo) 3/2 

Po->0 



-4/3 



dt 



t 3 Vt^T 



(83) 



(84) 



(85) 



(86) 



The integral is finite, so L(0) — 0, as in Figure 6 and Figure 8 



We now turn to the behaviour of strings with large po- For large p, the metric functions are 

e 4{*-».) = j _ ^ e -4 P /3 + 3 r _ a _ + 29 2 , e _ 8p/3 + (e _ 4p) (8?) 
c locr 
3iV f 



£ 2k = 2C g4p/3 



where 



This results in 



3 e 2$ ° 
2^72" 



f = {wf™^ + 32^ [ (2iVc " ^ f)2(1 ~ 8P) + 216N ^ £ ~ 4P/3 + ° ( 6 ~ 8P/3 ) 



1 2c 



e *oog2p/3 



, 3iV f 

1 

8c 



(89) 

(90) 
(91) 
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We again need the function tg/t'f, which for large p is 
t'f 2 V 64c^Vf 



[(2N C - N f ) 2 (8p+ 11) - 2497V f 2 ] e~ 4p/3 + O f^e" 8 "/ 3 ") 



Using ( 49 1 , for large p the separation therefore approaches 

T 



Nf. 



If we include the next term of the expansion ( 92 1 in the integral ( 35 ) , we get 



f 

L(p ) = Loo -2j 



dt 



tVt 2 



1 



where 



(92) 



(93) 



(94) 



(95) 



This means that L approaches L M from below. Together with the fact that the smooth solution 
has L(0) = 0, the result is behaviour which is qualitatively that of Figure 8 (c-f). 

For the behaviour of E, we need the function A'(T e g), defined in (40 1. In this case, for large 
p we simply have 



A'(T eff ) 



-e(p(T cS ))+0(e-^ 3 ) 

3(2iV c - iV f ) 2 ^3~logT (li 
8N f 



2c T 2 ff 



O 



1 

crt 



(96) 



This satisfies the condition (44 1, so the energy approaches a constant for large p. More precisely 

3(2N C -N f ) 2 /T 



( 43 1 results in 

E(p ) = £ 



where 



is a constant. 



47V £ 



2A(0) 



2c T cS (po) \T c ff(pa) 



dT r , 



(97) 



(98) 



The numerical calculation (Figure 12 1 confirms this, and reveals that a local maximum and 
minimum occur in L(p) for small JVf, while for large iVf we find that L(p) is always increasing, so 
that E(L) is smooth. In all the cases calculated E > for all L, and so the stable configuration 
is presumably the 'free' solution (251. 

Note that because a finite upper limit (in this case p\ = 20) is needed for the integrals, the 
numerical calculations cannot be trusted for large p. In particular, the numerical integration of 
(47 1 will always yield L(pi) — 0. The plots in Figure 12 have been terminated before L decreases 
significantly away from L^. 

By evaluating (23 1 numerically, we can determine the shape of the string. This is shown in 
|Figure 13) Of the strings shown all are unstable, except that with N{ = 5, po = 1, which falls 
within the region between the cusps in E{L) (Figure 12 1 and so is metastable. 

6.2.2 D3 brane 

The 't Hooft loop can be obtained by replacing the string in the above discussion with a D3 brane 
[271. with 



(p = 2n — (p, 



%l> = 71". 



(99) 
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(a) 



(b) 




10 



A) 




Figure 12: The results of the numerical calculation of (a) L(po) and (b) E(L) for a string in the 
'rotated' background, for various values of Nf. Here N c = 10 and hi = 25, resulting in values of c 
in the range 2.8 < c < 5.5, depending on Nf. 



(a) 



(b) 





Figure 13: The shapes of strings with varying po- Here the parameters are as in Figure 12 with 
(a) Nf = 5 and (b) Nf = 20. 
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As described inlsection 21 we therefore need to calculate 



where 



T cS (p) = T D3 e-*/ J d 2 X yJdetGfl, 



_(a)_ dO* dV 

U ab - 9ij 



(100) 



(101) 



' dX a dX b 

is the internal part of the induced metric on the D3 brane. If we choose the parametrisation 



{X 2 = 0, X 3 = tp}, we obtain from (69 1 



G™dX a dX b = /i 1/2 e 2 * 



e 2 "+-e 2 \a-l) 2 
4 



(d9 2 



V). 



Setting Td3 = 1 for convenience, this results in 



T eff = 47re q 



e 2ci +-e 29 {a-lf 
4 



(102) 



(103) 



Again using the asymptotic expansions in PQ , this is 



S7re*V 2 



2N{ 



logp 



1/4 



1 + 



log(-log p) 



2^ C e*~e 4p / 3 {l + - 
4c 



(2iV c -iV f )p + 



logp 
3A^ 



0. 



(104) 



' 4p/3 + 0(e- 2p ) 



As T c g(0) — we generically expect the smooth solution to have L(0) = 0, as in the case of the 
string, but to be sure it would again be necessary to evaluate the integral ( [82] ) . 

To obtain the large po behaviour, we again need the p-dependence of T c gg /T' cS f , and in this 
case we get an asymptotically constant separation 



37T r— 
Loo - -J-V N U 



(105) 



which is equal to half that obtained for the string. 



The results of the numerical calculations are shown in [Figure 14| The behaviour is qualitatively 



similar to that of the string (Figure 12 1, with the exception of the fact that in this case L approaches 



Loo from above. Although this appears to be an insignificant difference when viewed in terms of 
L(p ), the effect is that E decreases for large po. This means that there is a (small) range of L 
for which the smooth configuration is stable. 



6.2.3 Dl brane 

For a Dl brane the effective tension is 



T e & = e f, 



(106) 



where we set Tdi = 1. By (87 1, we see that for large p the tension approaches a constant multiple 



of that of the fundamental string, 



J eff ^ e 1 cS 



(107) 



The asymptotic value of L depends on the ratio T e g/T^ s and is therefore the same as in the case 
of the string, 



3tt 



Nf. 



(108) 
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Figure 14: The results of the numerical calculation of (a) L(po) and (b) E(L) for a D3 brane in 
the 'rotated' background. The parameters have the same values as in |Figure 12] Note that the 
endpoints of the curves were determined by the limitations of the numerical calculation. 



For p — > 0, we can write 



O 



( e 4(*-*-)) 



which results in 



Toff = 1 



N{ 



logp 



-3/2 



O 



log(- log p) 

(-l0gp)5/ 2 



(109) 



(110) 



As T c ff(0) ^ 0, the integrand in (47 1 is non-zero for p ^ 0, meaning that L(0) ^ 0. The contribution 
from the lower limit is of the form 

[dp (-logp)- 3/2 , (111) 
Jo 

which is finite. The separation is therefore finite and non-zero for po = 0, and the 'free' solution 



(25 1 does not exist. It is unfortunately not possible to determine L(0) analytically, because it 



would be necessary to evaluate the integral (47 1 over the whole range < p < oo 



The results of the numerical calculations are shown in |Figure 15] The behaviour for large po 
is confining, as expected when T e g ^ 0. 



(112) 



6.3 The 'unrotated' case 

In the case with f} — > 0, the metric is simply 

ds 2 = e* (dx 2 3 + dsg) , 

and / = e*. However, generically <!> — > constant for large p, and so the boundary condition ( p9| 
is not satisfied without /3 — > oo. This problem can be overcome by taking the limit c — > 0, which 
can be considered the flavoured generalisation of the solutions discussed in [2]. This case was 
discussed in |27j . from which we obtain the asymptotic expansion for large p, 



P = \2N C - Nf\ p 



1 



1 



2\2N C -N { \ p 



O 



(113) 



For small p the solution is unchanged from (|74|). Using (|72|) this results in 



/ 



A 



V2 



V|27V c -iVf|+iV f p- 1/4 e" [1 + Oip- 1 )] , 



(114) 
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Figure 15: The results of the numerical calculation of (a) L(po) and (b) E(L) for a Dl brane in 
the 'rotated' background. The parameters have the same values as in |Figure 12| 



where 



j i = \(\2N c -N { \+N { ) 2 \2N c -N f \ 



In this limit we still have / — > oo for p — !> oo, so (29 1 is satisfied 



(115) 





Figure 16: The behaviour of (a) P(p) and (b) Q(p) with c — > 0, showing the linear behaviour 
Q113||114| ) for large p. 



For the purposes of our discussion, the only difference between the 'rotated' and 'unrotated' 
backgrounds is 



,-2 = -2# 

•> unrotated ' 



f-2 _ -2$ _ -24> 
J rotated — c ' 



(116) 



For small p we have <& — > — oo, and so / is unchanged by the 'rotation'. The discussion of the 
previous section resulting in L(0) = for the string and D3 brane therefore also applies in the 
'unrotated' case. 



6.3.1 Fundamental string 



Using (1141, and the fact that for a fundamental string T c ff = /, we find that for large p, 



ff = % = ^V|2iV c -iV f |+iV f [1 + O(p-i)] 



(117) 
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so that the separation is asymptotically 



^=y/\2N c -N t \ + N t . 



(118) 



The numerical calculation again results in the expected modified 'soap film' behaviour (as in 



Figure 8 1, as shown in Figure 17 



(a) 



(b) 



10 

8 
6 
4 

2 [ 



0.5 



-0.5 



Po 

10 -1.0 



8 It 



Figure 17: The results of the numerical calculation of (a) L(po) and (b) E(L) for a string in 



the 'unrotated' background, with N c = 10 and N{ = 15. To obtain the linear behaviour (1131 
hi « 11.198, so as to give c = 0. 



6.3.2 D3 brane 

The internal space is unaffected by the 'rotation', so the D3 brane effective tension is still given 



by (103 1. However, we now have the limit c — > 0, so the large-p asymptotic expression (104) is no 

I O ( 'J . (119) 



longer valid. Instead, using the asymptotic solutions in 



T eff = 2wAe 9 ° \2N C - N t \ p 3/ V 



Together with (1141, this results in 



Lea — 



V2 



^/\2N C - N f \ + N f , 



(120) 



as for the fundamental string. 



6.3.3 Dl brane 

In the 'unrotated' background a Dl brane has constant tension 

T e s = T m e~^f — T D1 . 



(121) 



The minimum-energy solution to (21 1 is then simply a string which does not descend into the 



bulk, resulting in confining behaviour E = Tdi£. 



6.4 Discussion 



Qualitatively, the clearest result of the preceding analysis is the unusual behaviour seen in the 
flavoured solutions for small po. For both strings and D3 branes we find the 'soap film' behaviour 



shown in Figure 6 the primary result of which is that at large separations the only solution to the 



23 




Figure 18: The results of the numerical calculation of (a) L(po) and (b) E(L) for a D3 brane in 
the 'unrotated' background. The parameters have the same values as in |Figure 17| 



equation of motion (211 is the 'free' solution (251. This not particularly unusual; it is also seen 
for ex ample, in the case of the AdS's-Schwarzschild x S 5 background discussed in subsection 5.2 
In section 3 we related this behaviour to the fact that the effective tension vanishes for Rq = Rq 



It seems reasonable to suppose that this change relative to the unfiavoured case is related to the 
introduction of the IR singularity. For example, in the case of the fundamental string we have 
Tcs(po) ~ e*(po) for small po, and we see that this effect follows immediately from the fact that 
the dilaton diverges, Q(p) — > — oo, for small p. 

The modified behaviour for large po manifests itself in the fact that the separation goes to a 
non-zero constant for large po — > oo. This is also simplest when considered in the context of the 
fundamental string, so that T c g oc /. In that case the condition (38) for constant separation in 
the UV becomes 



9(r) 
f'(r) 



constant 



for r 



(122) 



for a generic radial coordinate r. 

condition becomes even simpler. From equations (90 1 and (114 1, we see that in the UV / is an 
exponential function of p, so that f'(p) 

g(p) 



f(p) 



In terms of the specific radial coordinate p used above the 
mations ( 90 1 and ( 114 1, we see tha 
f(p), and the relevant condition is 

for p -> oo. (123) 



constant 



It is clear that the constant-separation behaviour results from a precise cancellation between the 
functions / and g. This can be viewed as the fact that the 'exponential' coordinate p becomes 



identical (up to a constant) to the rescaled coordinate R introduced in subsection 2.2 

As taking p — > oo does not result in L(p ) — > 0, we find that the UV of the field theory 
(small separations) is no longer described by the large-/? region of the bulk theory. This is in 
contrast to the normal behaviour, and in particular that in the unfiavoured case, where increasing 
po corresponds to decreasing L. 

Instead, for separations less than some critical value the 'free' solution is stable, as is the case 
for large separations. In particular L = corresponds to a string reaching straight down to p = 0, 
which can be considered a degenerate case of both the 'free' and smooth solutions. Aside from 
this, the smooth solution describes at most only a small range of separations, and in many cases 
is unstable for all L. This describes non-interacting particles. 

However, it is not clear how much of this behaviour is physical. In particular, we might expect 
that resolving the IR singularity would result in L(Rq) ^ 0, presumably giving a confining IR 
as in the unfiavoured case (see the 'van der Waals' case, Figure 5 1. Unless we also then have 
L — > for R — > Roo we would gcncrically expect a minimum separation, which would be difficult 
to understand physically. 



24 



Acknowledgments 

I would like to thank C. Nunez and G. Silva for valuable discussions. This work was supported 
by an STFC studentship. 



25 



References 



[i] 

[2] 
[3] 
[4] 
[5] 

m 

[7] 
[8] 

[9] 
[10] 

[11] 
[12] 

[13] 
[14] 
[15] 

[16] 
[17] 

[18] 

[19] 

[20] 
[21] 
[22] 
[23] 
[24] 
[25] 

[26] 
[27] 



J. Gaillard, D. Martelli, C. Nunez, and I. Papadimitriou (2010), |arXiv :1004.4638 
K. G. Wilson, Phys. Rev. D 10, 2445 (1974). 



J. M. Maldacena, Adv. Theor. Math. Phys. 2, 231 (1998), |arXiv:he p-th/9711200 
J. Sonnenschein (1999), |arXiv:hep-th/0003032 



S.-J. Rey and J.-T. Yee, Eur. Phys. J. C22, 379 (2001), |arXTv:hep-th/9 803001 



J. M. Maldacena, Phys. Rev. Lett. 80, 4859 (1998), arXiv:hep-th/9803002 



C. Nunez, M. Piai, and A. Rago, Phys. Rev. D81, 086001 (2010), arXiv:0909.0748 



N. Drukker, D. J. Gross, and H. Ooguri, Phys. Rev. D60, 125006 (1999), [arXivihep- 
|th/9904191| 



A. Brandhuber and K. Sfetsos, Adv. Theor. Math. Phys. 3, 851 (1999), arXiv:hep-th/9906201 



S. D. Avr amis, K. Sfetsos, and K. Siampos, Nucl. Phys. B769, 44 (2007), |arXiv:hep- 
|th/0612139| 

C. Bachas, Phys. Rev. D 33, 2723 (1986). 



Y. Kinar, E. Schreiber, and J. Sonnenschein, Nucl. Phys. B566, 103 (2000), arXiv:hep- 



th/9811192 



I. R. Klebanov and M. J. Strassler, JHEP 08, 052 (2000), arXiv:hep-th/0007191 



J. M. Maldacena and C. Nunez, Phys. Rev. Lett. 86, 588 (2001), [arXivlhep -th/0008001 

F. Bigazzi, A. L. Cotrone, C. Nunez, and A. Paredes, Phys. Rev. D78, 114012 (2008), 
arXiv: 08007411 



F. Bigazzi, A. L. Cotrone, and A. Paredes, JHEP 09, 048 (2008), arXiv:0807.0298 



F. Bigazzi, A. L. Cotrone, and A. Paredes, Nucl. Phys. Proc. Suppl. 192-193, 134 (2009), 
larXiv:0810.40l8l 

F. Bigazzi, A. L. Cotrone, A. Paredes, and A. V. Ramallo, JHEP 03, 153 (2009), 
larXiv:0812.3399l 

C. Nunez, I. Papadimitriou, and M. Piai, Int. J. Mod. Phys. A25, 2837 (2010), 
larXiv:0812T3655l 

S. D. Avramis, K. Sfetsos, and K. Siampos, Nucl. Phys. B793, 1 (2008), [arX"iv:0706.2655| 



E. Witten, Adv. Theor. Math. Phys. 2, 505 (1998), arXiv:hep-th/9803131 



R. E. Arias and G. A. Silva, JHEP 01, 023 (2010), arXiv:0911.0662 



J. Maldacena and D. Martelli, JHEP 01, 104 (2010), arXiv:0906.0591 



R. Casero, C. Nunez, and A. Paredes, Phys. Rev. D73, 086005 (2006), arXiv:hep-th/0602027 



A. Butti, M. Grana, R. Minasian, M. Petrini, and A. Zaffaroni, JHEP 03, 069 (2005), 
|arXiv:hep-th/0412187| 



R. Casero, C. Nunez, and A. Paredes, Phys. Rev. D77, 046003 (2008), arXiv:0709.3421 



C. Hoyos-Badajoz, C. Nunez, and I. Papadimitriou, Phys. Rev. D78, 086005 (2008), 
arXiv:0807.3039. 



26 



